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ABSTRACT 



The direction of flow along the centerline of several sections 
of the inlet duct of a centrifugal pump was determined using a hot wire 
anemometer. In contrast to other tests of this kind, the pump was 
operated with air to facilitate the measurements. The present investi- 
gation is part of a program at the U, S. Naval Postgraduate School to 
evaluate the performance of liquid pumps with air in order to apply 
aerodynamic measuring methods and equipment to the task. Data were 
taken with the pump operating at different points on the characteristic 
curve. The flow direction data are presented as a fxinction of position 
in the intake duct for various operating conditions. 

Prerotation was observed, the character of which was in agreement 
with theoretical predictions in that prerotation was in the same direc- 
tion as the impeller rotation at the low flow rates and in the opposite 
direction at the high flow rates. Backflow near the periphery of the 
duct was observed at the low flow rates. 
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AN INVESTIGATION OF THE PREROTATION CHARACTERISTICS 



OF A FLUID IN THE INTAKE DUCT OF A CENTRIFUGAL PUMP 

1. Introduction. 

It has been observed for many years that rotation of fluid in the 
inlet ducts ahead of an impeller exists. This rotation may be in the 
same direction as, or opposite to, the direction of rotation of the 
impeller. The effect is more pronounced and easily observable at off 
design conditions but may even occur at the design point. It is essen- 
tial to know the fluid angles at the impeller inlet in order to properly 
design the inlet blade angles. If the blade angle is very different 
from the relative flow angles, local high fluid velocities may produce 
the phenomenon of cavitation. While some cavitation may be tolerated 
to achieve higher flow rates, it is generally undesirable because it 
is destructive and causes inefficient and unstable operation of the 
pump. 

The cause of the phenomenon o prerotation is somewhat obscure and 
little quantitative data is available to pump designers. Although con- 
ventional hydraulic pumps have been examined for prerotation, there do 
not exist prerotation data in the literature for the very small blade 
eingles encountered in rocket fuel and oxidizer pumps. Ebcperiments that 
have been done generally have used liquid pumps with water as the working 
fluid. It will be seen that the same results may be obtained using a 
liquid pump but with air as the working fluid. The latter wotild require 
much less power to drive it, less elaborate experimental equipment, and 
would provide easier ways to investigate the physical aspects of the process. 

2. The Pump. 

The pump used for this experiment was a liquid fuel pump for the 
Titan I missile. It is a centrifugal pump seven inches in diameter at 
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the eye and 11 inches at the discharge of the impeller. It has a diffuser 
with four vanes, an anniolar collector and a three inch diameter discharge 
flange. There is a static pressure tap at the discharge flange where 
the pressure, P 2 , was taken. The impeller has two blades at the eye 
with additional blades starting at larger radii for a total of eight 
blades at the discharge. Figs. 1, 2, 3. The pump was driven by an air 
turbine , 

3, The Inlet Duct. 

The inlet duct consisted of a seven inch diameter, four inch long 
Plexiglas spacer, a seven inch diameter, 24 inch long steel tube and 
a built up wooden bell mouth. Through the plexiglas, qualitative 
observations of the flow were made using smoke and cotton tufts. Mounted 
on the steel tube were four brackets with a hole through the pipe 'tall, 
through which a probe coixLd be inserted to make a survey of the flow 
along the centerline of the pipe, top to bottom. These survey points 
were six inches apart, with the nearest one to the impeller being seven 
inches from the inlet flange. This point is referred to as position 
number four. 

4. The Discharge Pipe. 

The discharge pipe was a three inch pipe, ten feet long. It had 
a bundle of one-half inch copper tubing six inches long at the pump end 
to act as a flow straightener. Five feet from the pump was a flow meter 
consisting of a total pressure Kiel probe with a thermocouple mounted at 
the centerline of the pipe and two static ports on opposite sides of the 
pipe. For calibration of the flow meter see Appendix C, At the discharge 
end of the pipe was a gate valve for controlling the flow rate. 
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5. The Hot Wire Anemometer. 

A Plow Corporation, Model HWB, hot wire anemometer was used to 
determine flow direction using a .00035 inch diameter tvingsten wire. 

The probe was motmted in a protractor which in turn was mounted on a 
slide. The depth of the hot wire could be determined by a scale on 
the slide and the directional orientation of the hot wire could be read 
to 0.2 degrees using a vernier scale on the protractor. The centerline 
of the inlet duct at position number two at the highest flow rate was 
used as a calibration point for directionally orienting the hot wire. 

6. Other Equipment. 

An electronic counter was used to determine pump speed. A bank of 
vertical tube manometers was used to measure pressures in cm. of water. 
The differential pressure was measured with a water micromanometer. 
Discharge total temperature was determined with a Brown thermocouple 
instrument. A four-foot cubic plywood box was used as an inlet plenum 
to act as a buffer between the inlet duct and ambient turbulence. The 
inlet to the plenum was a two-foot square hole covered with four layers 
of filter cloth. 

7. Speed of the Pump. 

The design speed of the pump is 14,000 RPM. It was originally 
planned that several speeds in the neighborhood of the design speed 
would be used for this experiment. Upon running the pump at speeds of 
10,000 to 12,000 RPM, large amplitude vibrations were observed. A 
dynamic balance of the air turbine failed to lessen the vibrations so 
a speed of 8,800 RPM was settled upon for the tests. The vibrations 
are evidently due to the ball bearings ^ose axial loading is not as 
great as design loading because air is the fluid being pumped. 
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In order to avoid cavitation the eudal velocity ahead of the pxamp 
is kept small. The reduction in RPM produced a further decrease in 
velocity which caused difficulty in measuring the dynamic pressure in 
the duct. 

8, Measuring the Flow Angle. 

The first attempt at measiiring the flow direction in the inlet duct 
\ms made utilizing a United Sensor airflow probe. This probe had five 
pressure ports. When the probe was properly oriented in the air stream, 
one would measure the total pressure, two would measure the static 
pressure, and two others not used in this application could be used to 
measure components of the flow in the direction of the probe axis. 

After calibration in an air stream of known direction, the probe could 
be used to determine the direction of airflow in the inlet duct. Not 
only the direction, but also the velocity of the flow could be determined 
at any point in the duct tested, and from this the flow rate could have 
been calculated. 

The very low flow rates made the pressvire measurements to determine 
flow angles very critical. The maximum flow rate which was measured 
later with a flow meter was .46 Ib/sec. This flow rate, assuming a 
uniform velocity distribution in the inlet duct, irould yield a dynamic 
head of about 0.1 inches of water. Consider now the ideal flow about 
a cylinder where 

Ps - Ft = qo(l-4sin2e) » O.ld-Asin^e) 

then dPs = -0.8 sinGcosQ 

d0 

for e *= 45° 

dPs - -0.4 in.H20/radian s -0.007 in.H20/degree 
de 
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Hence, if the probe were one degree away from facing directly into the 
stream, the pressure difference between the two static ports x/ould be 
about ,015 inches of water. An inclined manometer was available which 
could be read to ,01 inches. However, since the sensitivity goes down 
even further as the flow rate is reduced, it was felt that the accuracy 
of results would be rather poor with this method. 

The hot wire anemometer offered a possible solution. As pointed 
out in Appendix D, the directional sensitivity of the hot wire is great- 
est when the wire axis is aligned with the direction of flow. It was 
felt that the supporting pin ahead of the wire might cause turbulence 
and shield the flow from the full effect of the cooling gas stream. To 
avoid this, the leading end support pin was filed off so that the wire, 
when soldered across the tips of the pins, made an angle of about 65*^ 
with the probe axis. No apparent reduction in the galvanometer needle 
fluctuations was noted. Since the directional sensitivity of the probe 
is reduced by this technique, the standard probe was used. 

The hot wire anemometer, being sensitive enough to detect very 
small velocities, is also sensitive to perturbation velocities. There- 
fore, it became necessary to build the intake plenum chamber, described 
in Sec. 6, to prevent ambient turbulence from giving an indication of a 
change in flow direction. Even with the plenum chamber there was always 
some fluctuation of the gsO-vanometer needle used to determine the resis- 
tances. With fluctuations of the needle it was difficult to find a 
maximum, so the scheme actually used was to record the two angles, 

15-20 degrees on either side of the angle for maximum deflection, which 
produced the same galvanometer reading. The angle halfway between the 
two recorded was then used as the flow angle. In this technique there 
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can be no shielding by the upstream supporting pin because it is not 
directly in front of the v/ireo 
9. Existing Theories. 

T'lhlle the exact mechanism of prerotation is not clear, the phenom- 
enon has been observed by many investigators for years. £61.. Several 
attempts have been made to explain the reasons for prerotation. At 
design conditions the fluid is said to acquire prerotation by taking the 
path of least resistance or least inlet impact. For flow rates lower 
than design the flow acquires prerotation which is the sane direction 
as the impeller rotation, a positive prerotation, and at flow rates 
higher than design the prerotation is negative. [5j . Backflow in the 
outer portion of the impeller channel at low flow rates, and high net 
heads are believed to be the cause of positive prerotation. Backflow 
has been observed in the inlet duct and in the diffuser. [4] . 

The impeller discharge conditions seem to have a great deal of 
influence on prerotation. Experiments at California Institute of 
Technology using a vaneless, symmetrical diffuser showed practically 
no evidence of prerotation over the entire range of flow rates. £2]. 

This pump did show a positive component of tangential velocity at the 
discharge of the impeller when operating at zero pressure coefficient. 
This positive angular momentum at the discharge would seem to indicate 
a positive rotation at the inlet, but no such rotation was observable. 

At zero pressure coefficient the accepted theory would predict a negative 
prerotation. Shepherd believes that the positive angular momentum at the 
discharge is due to losses in the impeller. 

The pump studied by Peck had a vaneless diffuser, but had the typi- 
cal eccentric discharge volute. [4] » This pump exhibited considerable 
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prerotation and backflow at low flow rates „ It showed no negative 
prerotation at flow rates above the design point o 

Existing theories, therefore, seem to be of little practical use 
to the ptimp designer. The discharge conditions evidently are a contrib- 
uting influence on prerotation „ 

10. Results. 

The results of these tests are tabulated in Tables 1 through 22 
and are presented in graphical form in Figs. 8 through 17. 

11. Discussion. 

Fig. 8 is the experimentsilly determined characteristic curve of 
the pump. The pressure coefficient used here is non-dimensional and 
is defined as follows; 



where, P 2 is a derived density using the static pressure at the outlet 
flange and the total temperature at the flow meter. The flow coefficient 
is also non-dimensional and is defined as follows; 



\2> being measured directly, is a derived quantity using the 
measured flow rate, w, the density, and outlet area of the impeller. 

The outlet area does not consider any restriction by the blades and is 
merely TTD 2 b„ 



The theoretical pressure coefficients as a function of the flow 
coefficient were computed using the slip factors of Stodola and Busemaan. 
These are plotted in Fig, 6 for reference. The Stodola slip factor, 
was calculated using his equation; £5j 





\ 2/»2 



V = ''//■amDjb 



iUg = 1 - TTC0^2B 



V^2 (theoretical) 



n 
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Busemann's slip factor, /ig, for eight blades and an exit blade angle 
of -73.75°, Fig. 3, was obtained from a graph in Ref. 5. The theoretical 
pressure coefficient then may be expressed ass 

"VaSjB =>“S,B ^ 

The efficiency may be expressed as the ratio of actual to theoret- 
ical pressure coefficients. The efficiency was calculated and plotted 
in Fig. 7. The pressure coefficients and efficiency can be expected to 
increase somewhat with the liquid fuel as a working fluid due to increased 
Reynolds number. 

The prerotation observed in these tests is generally what xiould be 
predicted by accepted theories. Positive prerotation and backflow in 
the outer portions of the inlet duct were observed at low flow rates. 
Positive prerotation is indicated in the figures as an angle greater 
than zero for depths less than three and one hadf inches and an angle 
less than 360 degrees for greater depths. 

There is considerable scatter in the data with respect to the flow 
angle. Several repeatability tests were made during the course of the 
experiment which indicated that the data should be accurate to within 
one degree. These tests were performed in connection with the direc- 
tional calibration of the hot wire. Each time the wire had to be 
replaced it was necessary to establish a new angular tare. Since this 
was a one time measurement which effected all the readings taken x^ith 
that x/lre, it was considered worthwhile to use an average of eight to 
ten readings rather than just one. No readings were more than one degree 
different from the average. However, these tests were made on the 
centerline of the inlet duct at position number two with the pump 
operating at maximum flow coefficient. Directional sensitivity of the 
hot wire is proportional to the velocity so that the accuracy would be 
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poorer at the lower flow rates. Velocity perturbations in any direction 
would give an indication of a change in flow direction to the observer. 

It can reasonably be expected that these perturbations would be 
oscillatory in character, however. 

Increasing the speed of the pump to 14,000 RPM wo\xld considerably 
increase the accuracy of the data. Directional sensitivity is propor- 
tional to the velocity with the hot vlre anemometer and proportional 
to the square of the velocity if the pitot-static tube is used. Velocity 
is proportional to the RPM, Perturbation velocities would probably not 
increase as much so the fluctuations of the needle would be reduced also. 

In spite of the scatter, it is felt that some negative prerotation 
was present at the high flow rates. In Fig, 17, an interpretation has 
been made which shows a maximum of two degrees of negative prerotation 
at the maximum flow rate. 

At the apparent design point there is still considerable positive 
prerotation and backflow. The point of highest efficiency is at a flow 
coefficient of .03, Fig. 12, showing data taken at a flow coefficient 
of .048, shows that positive prerotation and backflow were observed at 
position number four. 

Velocity measurements were not taken so the exact incidence angles 
cannot be calculated. However, by assuming that the magnitude of the 
velocity is a constant across the section for each flow rate, an estimate 
of the effect of prerotation on the incidence angles may be made. This 
was done for the two extreme flow rates. The follo\d.ng analysis was 
made to determine the velocities in the inlet duct for the case with 
prerotation and for the case of axial flow with uniform velocity 
distribution for comparison. 
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T = = 2TTv/r cos c<]^ dr 



- 2TTV^rj cosc<2^jAr^ 

j 

for the axial flou case 

= TTY{t/ - ri^) 



The results were; 

i >2 = *026 



V(wlth rotation) = ^'Vsec 

axial flow) ~ 7«1 ft/sec 

<j >2 = «097 ^(with rotation) = ft/sec 

axial flow) = 26.4 ft/sec 

These results were used in plotting the velocity triangles in Figs. 19 
through 22. 

It can be seen that, except for the very large amounts of prero- 
tation at the low flow rate, velocity of the flow has a greater influence 
on the incidence angles than the flow angle itself. Of course, these 
velocities must be greater with prerotation than without for any given 
flow rate. Therefore, the tendency to cavitate would be greater with 
prerotation because the increased velocities produce reduced static 
pressures. 

12 . Conclusions . 

These tests show essentially the same trands as other tests have 
shown where water was used as the working fluid. 

The very small velocities involved may be handled very well with 
the hot wire anemometer. This instrument may be used to measure the 
velocity as well as the flow direction. Appendix D. 

Testing liquid pumps with air has many advantages. The power 
requirements for driving the pump are less. Return and supply lines 
with their accompanying piimps and coolers are not required. Measuring 
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equipment is not subjected to high dynamic loads. Test setups may be 
kept simple and aerodynamic measuring devices can be used. Packing is 
not required aromd the probe to prevent leakage and the unused holes 
in the inlet duct may be sealed with a piece of drafting tape, 

13 . Recommendations , 

It is recommended thats 

a. Efforts be made to eliminate vibrations so that tests 
may be made at higher speeds. Preloading the bearings 
and a dynamic balance of the entire assembly should reduce 
vibrations. Installation in one of the new test cells in 
a more permanent fashion would eliminate feedback of 
vibrations from the supply line and the discharge pipe. 
Directional sensitivity of the hot wire is proportional 

to the velocity so this would improve the accuracy of the 
data, 

b. Velocity measurements as well as direction measurements 
be taken using the hot wire anemometer. Direction should 
be determined first and then velocity should be measured 
with the wire axis perpendicular to the air stream, 

c. A more sophisticated probe which employs two hot wires 
on the same probe be purchased or fabricated to try to 
obtain more accurate data, 

d. The diffuser vanes be removed to determine the effect of 
that section of the pump on the prerotation, 
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TABLE 1 



RUN NO. 1 



FLOW COEFFICIENT 



PREROTATION TEST 



12 MARCH 1964 

POSITION NO. 1 j 

. ' ' « ' 

* .07594' PRESSURE COEFFICIENT = .27850 



UEPTH 


ANGLE 


.10 


355.40 


.50 


2.60 


1.00 


1 .60 


1.50 


357.10 


2.00 


357.80 


2.50 


358.70 


3.00' 


359.20 


3.50 


.20 


4.00 


359.60 


4.50 


1 .40 


5.00 


1.80 


5.50 


.30 


6.00 


3.80 


6.50 


1 .30 


6.90 


356.50 
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table 2 



PREROTATION TEST ' , 

RUN NO. 2 • 13 MARCH 1964 . 

POSITION NO. 1 



FLOW COEFFICIENT = .09598 

CEPTH 

.10 

.50 

1.00 

1.50 

2.00 

2.50 

3.00 

3.50 

4.00 

4.50 

5.00 

5.50 

6.00 

6.50 
6.90 



PRESSURE COEFFICIENT = .08463 

I 

) 

ANGLE . ■ ■ , 

.90 

1.20 

2.00 

355.60 
357.00 

359.60 

2.70 

• I 

.00 

3.70 . *’ 

359.40 

3 . 80 . 

1.10 

5.50 

355.40 
355.70 
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TABLE 3 



RUN NO. 3 



FLOW COEFFICIENT 



PREROTATION TEST, 



13 MARCH 



, POSITION NO. 4 

= .09590 PRESSURE COEFFICIENT * 



oepth 


ANGLE 


.10 


1.90 


.50 


359.00 


1.00 


359.10 


1.50 


353.95 


2.00 


358.50 


2.50 


357.30 


3.00 


359.10 


3.50 


.95 


4.00 


.50 


4.50 


2.15 


5.00 


1.50 


5.50 


1 .10 


6.00 


.20 


6.50 


358.70 


6.90 


359.20 



1964 



08522 
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TABLE k 



PREROTATION TEST 

RUN NO. 4 13 MARCH 1964 

POSITION NO. 4 



FLOW COEFFICIENT = .06336 

DEPTH 

.10 

.50 

• 1.00 
1 .50 
2.00 

2.50 

3.00 

3.50 

4.00 

4.50 

5.00 

5.50 

6 . 00 

6.50 
6.90 



PRESSURE COEFFICIENT = .37310 

ANGLE 

1.60 

359.70 

.45 

352.80 

I 

357.00 

357.20 

358.95 

1.50 

.35 

359.40 

358.75 

.30 

) 

•.15 

359.90 

356.55 
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TABLE 5 



RUN NO. 5 



FLOW COEFFICIENT 



PREROTATION TEST • 

13 MARCH 1964 

POSITION NO. 4 

= .02145’ PRESSURE COEFFICIENT = .56532 



DEPTH 


ANGLE 


. 1,0 


102.50 


.50 


97.00 


1.00 


84.60 


1.50 


67.20 


2.00 


45.90 


2.50 


27.95 


3.00 


12.20 


3.50 


1.30 


4.00 


353.30 


4.50 


340.45 


5.00 


322.40 


5.50 


299.40 


6.00 


280.20 


6.50 


264.95 


6.90 


256.10 
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TABLE 6 



RUN NO. 6 



FLOW COEFFICIENT 



PREROTATION TEST 

» 

17 MARCH 1964 

POSITION NO. 4 

= .0480'8 PRESSURE COEFFICIENT = .49974 



CEPTH 


ANGLE 


.10 


101.45 


.50 


95.00 


1.00' 


78.25 


1.50 


55.70 


2.00 


21.50 


2.50 


9.50 


3.00 


2.80 


3.50 


360.00 


4.00 


358.00 


4.50 


354.10 


5.00 


341.40 


5.50 


321.00 


6.00 


289.70 


6.50 


267.75 


6.90 


252.30 
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TABLE 1 . : 

PREBOTATION TEST 

RUN NO. 7 17 MARCH 1964 

I ' 

\ ■ , 

POSITION NO. 4 

FLOW COEFFICIENT = .06158' PRESSURE COEFFICIENT = .40929 



DEPTH 

1 


ANGLE 


. 10 


357.50 


.50 


357.60 


1.00 


354.90 


1.50 


355.60 


2.00 


357.30 


2.50 


357.30 


3.00 


356.50 


3.50 


.05 


4.00 


359.65 


4.50 


354.00 


5.00 


359.35 


5.50 


358. UO 


6.00 


359.70, 


6.50 


355.40 


6.90 


351.90 



21 




I 

I 

I 



I 

I 



1 



TABLE 



8 



RUN NO. 8 



FLOW COEFFICIENT 



PREROTATION TEST, 

18 MARCH 196U 

POSITION NO. U 

» 

= .07989 PRESSURE COEFFICIENT = .27829 



CEPTH 


ANGLE 


.10 


.15 


.50 


;40 


1.00 


359.50 


1.50 


. 10 


2.00 


358.85 


2.50 ■ 


359.90 


3.00 


.45 


5.50 


.40 


4.00 


.65 


4.50 


358.85 


5.00 


359.50 


5.50 


359.95 


6.00 


1.15 


6.50 


2.05 


6.90 


1.55 
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TA3LE 9 



PREROTATION TEST 

RUN NO. 9 18 MARCH 1964- 

POSITION NO. 4 , 

FLOW COEFFICIENT * .06496 PRESSURE COEFFICIENT = .40275 



DEPTH 


ANGLE 


: 10 


3.45 


.50 


2.40 


1.00 


.90 


1.50 


358.80 


2.00 


359.60 


2.50 


359.35 


3.00 


358.90 


3.50 


359.70 


4.00 


1.00, 


4.50 


.50 


5.00 


358.60 


5.50 


1.45 


6.00 


.75 


*6.50 


2.80 


6.90 


.90 
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1 

I 



14 




TABLE 10 



PREROTATION TEST ' . ' 

RUN NO. 10 ' 18 MARCH 196U 

POSITI ON NO. 3 i , 

FLOW COEFFICIENT = .065L3 PRESSURE COEFFICIENT = .1<0C59 



DEPTH 


ANGLE 


.10 


2.60 


.50 


3.20 


1.00 


.70 


1.50 


358.55 


2.00 


359.40 


2.50 


357.05 


3.00 


359.15 


3.50 


359.55 


U.OO 


359.10 


4.50 


1.95 


5.00 


.80 


5.50 


.50 


6.00 


359.60 


6.50 


1 .05 


6.90 


.50 



2h 



k 

r 

l 

I 

I 

I 

i 

j 









TABLE n 



PREROTATIGN TEST 

' “A 

RUN NO. n 18 MARCH 1964,. 

POSITION NO. 3 

FLOW COEFFICIENT = .0970?' PRESSURE COEFFICIENT = .08362 



DEPTH 


ANGLE 


. 10 


1.10 


.50 


2.60 


1 .00 


.50 


1.50 


358.25 


2.00 


359.10 


,2.50 


359.00 


3.00 


359.60 


3.50 


.20 


4.00 


1.60 


4.50 


1.30 


5.00 


1 .35 


5.50 


‘359.30 


6.00 


.45 


6.50 


359.30 


6.90 


.75 
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TABLE 12 



, PREROTATION TEST 

4 

RUN NO. 12 19 MARCH 1964 

f 

POSITION NO. 1 



FLOW COEFFICIENT = .08693 

DEPTH 

.10 

.50 

1.00 

1.50 

2.00 

2.50 

3.00 

3.50 
4. CO 

• 4.50 
’V 5.00 

5.50 

6.00 

6.50 
6.90 



PRESSURE COEFFICIENT = .20918 
ANGLE ' • 

4.10 

1 .00 

! 

I 

.50 

357.40 

2.10 
357.00 
358.20 

359.50 

359.50 

.20 

.10 

359.50 ' • 

358.10 

359.05 

358.50 
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TABLE 13 



RUN NO. 13 



FLOW COEFFICIENT 



PREROTATION TEST 



24 MARCH 1964 



POSITION. NO. 1 

= .03853 PRESSURE COEFFICIENT = .56494 



DEPTH 


ANGLE 


, .10 


5.70 


.50 


5.50 


1.00 


3.00 


1 .'50 


2.20 


2.00 


1.20 


2.50 


355.80 

^ ’ 


3.00 


352. ,50 


3.50 


353.30 


.4.00 


357.80 


4.50 


357.90 


’5.00 


359.80 


5.50 


.50 


6.00 


.70 


6.50 


359.60 


6.90 


355.30 
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T AB L E 1 U 



RUN NO. 14 



FLOW COEFFICIENT 



PREROTATION TEST 

31 MARCH 1964 

POSITION NO. 2, 

= .09675 PRESSURE COEFFICIENT = .08516 



DEPTH 


ANGLE 


• .10 


2.40 


.50 


4.50 


1.00 


2.20 


1.50 


358.00 


2.00 


356.30 


2.50 


358.60 


3.00 . 


359.70 


3.50 


.00 


4.00 


1.35 


4.50 


|'.20 


5.00 


3.00 


5.50 


1 .80 


6.00 


359.40 


6.50 


358.70 


6.90 


357.00 
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TABLE 15 



PREROTATION TEST 

RUN NO. 15 . 31 MARCH 1964 

tl 

POSITION NO. 2 

FLOW COEFFICIENT = .0792?’ PRESSURE COEFFICIENT = .28404 



DEPTH 


ANGLE 


.10 


358.70 


. 50 


358.90 


1.00 


357.70 


1.50 


.15 


2.00 


358.45' 


2.50 


358.60 


3.00 


.10 


3.50 


358.60 


4.00 


.80 


4.50 


1.50 


5.00 


1.50 


5.50 


1.30 


6.00 


359.20 


6.50 


358.70’ 


6.90 


.45 
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TABLE 16 



RUN NO. 16 



FLOW COEFFICIENT 



PREROTATION TEST 

1 APRIL 

POSITION NO.. 2 

= .06435 PRESSURE COEFFICIENT = 



DEPTH 


ANGLE 


.10 


3.00 


.50 


1 .90 


1.00 


359.35 


1.50 


358.30 


2.00 


358.40 


2.50 


358.90 


3.00 


.70 


3.50 


1 .05 


4.00 


.40 


4.50 


.75 


5.00 


360.00 


5.50 


1 .50 


6.00 


359.50 


6.50 


.90 


6.90 


1 .30 



1964 



40481 
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TABLE 17 



RUN NO. 17 



FLOW COEFFICIENT 



PREROTATION TEST 

1 APRIL 

. POSITION NO. 2 

= .04756 PRESSURE COEFFICIENT » 



DEPTH 


ANGLE 


.10 


4.90 


.50 


359.60 


1.00 


.60 


1.50 


358.10 


2.00 


357.70 


2.50 


357.90 


3.00 


360.00 


3.50 


359.80 


4.00 


.30 


4.50 


.40 


5.00 


1.35 


5.50 


.30 


6.00 


.•70, 


6.50 


359.40 


6.90 


3.30 



1964 



52388 



31 



s 



\ 

I 

i 




TABLE 18 



RUN NO. 18 



FLOW COEFFICIENT 



PREROTATION TEST 

1 APRIL 

POSITION NO. 2 

= .02765 PRESSURE COEFFICIENT = 



DEPTH 


ANGLE 


.10 


104.30 


.50 


100.30 


1 .00 


87.10 


1.50 


66.20 


2.00 


43.10 


2.50 


22.60 


3.00 


1 1.50 


3.50 


2.20 


4.00 


357.40 


4.50 


346.50 


5.00 


330.50 


5.50 


306.30 


6.00 


282.60 


6.50 


264. 10 


6.90 


256.60 



1964 , 



60983 
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TABLE 19 



PREROTATION TEST 

I 

RUN NO. 19 1 APRIL 1964 

POSITION NO. 3 



FLOW COEFFICIENT = .02663 

DEPTH 

.10 
.50 
1.00 
1 .50 
2.00 

2.50 

3.00 

3.50 

4.00 

4.50 
,5.00 

5.50 

6.00 

6.50 
6.90 



PRESSURE COEFFICIENT = .60978 

ANGLE 

103.80 

96.80 

87.30 

74.20 

56.20 

I ‘ 

’ 35.90 

18.10 
8.70 

352.60 

334.90 • ; 

•ir t - , r 

il3^60 ( • 

295.85 

277.60 
263.40 

256.60 
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TABLE 20 



RUN NO. 20 



FLOW COEFFICIENT 



PREROTATION TEST 

I 

1 APRIL 

f 

POSITION NO. 3 

= .07792 PRESSURE COEFFICIENT = 



DEPTH 




ANGLE 

t 


. 10 




2.40 


.50 




2.10 


1.00 




359.50 


1.50 




359.00 


2.00 


1 


358.35 


2.50 




.20 


3.00 




.60 


3.50 




.30 


4.00 




1 .30 


4.50 




.60 


5.00 




359.70 


5.50 


' 


'.80 


6.00 




359.80 


6.50 




360.00 


6.90 




359.60 



, * 
1964 



29394 






TABLE 21 



RUN NO. 21 



FLOW COEFFICIENT 



PREROTATION TEST 

2 APRIL 

POSITION NO. 3 

= .08694 PRESSURE COEFFICIF.NT = 



DEPTH 


ANGLE 


• . 10 


.10 


.50 


359.80 


1.00 


357.50 


1.50 


358.20 


2.00 


359.00 


2.50 


.60 


3.00 


360.00 


3.50 


.70 


4.00 


1 .00 


4.50 


1 .60 


5.00 


1.10 


5.50 


360.00 


6. 00 


359.50 


6.50 


359.70 


6.90 


.20 



1964 



19722 
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TA3LE 22 



PREROTATION TEST 



RUN NO. 22 2 APRIL 

POSITION NO. 4 

FLOW COEFFICIENT = .08687 PRESSURE COEFFICIENT = 



DEPTH 


ANGLE 


. 10 


.90 


.50 


356.00 


1.00 


.50 


1.50 


358.60 


2.00 


359.00 


2.50 


359.10 


3.00 


359.50 


3.50 


1 .10 


4.00 

% 


,.70 


4.50 


.10 


5.00 


359.60 


5.50 


359.10 


6.00 


.60 


6.50 • 


357.90 


6.90 


2.10 



1964 



19671 
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APPENDIX B 
FIGURES 
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/IG. 1, Photographs of Inpeller 
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FIG. li. Photograph of the Inlet Duct 
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PIG. 5* Photograph of the Flow Meter 





Fig. 6. Photographs of 
3ntire Test Setup 
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APPENDIX C 

CAUBRAnON OF THE FLOW METER 
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APPENDIX C 



CALIBRATION OF THE FLOW METER 



Since flow rates were very small without obstnictions in the 
discharge line, it was felt that some method of determining flow rates 
with a minimum obstruction would best suit the purposes of this experi- 
mento The only obstruction introduced into the flow in the flow meter 
devised was a Kiel probe. For purposes of calibration the discharge 
pipe was connected to an air supply line in which was installed a 
standard sharp-edged orifice type flow meter, so that all the air which 
passed through the standard meter also passed through the discharge pipe. 
The static pressure, P|, the difference between the centerline total 
pressure and the static pressure,^', and the centerline total tem- 
perature, T 2 , were measured in the discharge pipe. Fig, 5. The 
following two equations were to be used for determining flow rate, 

(1) PgV^ . fAP' 

2gRT2 

(2) w = P 3 A7 

RT2 



combining, 

(3) w = A( 2 AP‘lgPg/RT 2 )^ 

Everything can be measured except the factor 'f . The object of the 
calibration then is to determine T . This factor was not a constant so 
following the usual practice in flow meters, was plotted as a function 
of Reynolds number. Fig. 18 . In this graph a modified Reynolds 
number, Rex, was used, where, 

Rex r Re/3600 = WP 

3600/1 
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APPENDIX D 

APPLICATIONS OF THE HOT WIRE ANEMOMETER 
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APPENDIX D 



APPLICATIONS OF THE HOT I^RE ANEMOMETER 

The hot wire anemometer operates on the principle that the resist- 
ance of most metal wire is a linear function of the temperature: 

R = Rail f (T - To)] 

where, Rq is the resistance at some reference temperature, T^; and (Xq 
is the ' temperature coefficient of resistivity at that temperature. 

It has also been found that, if a very fine wire is heated by an 
electric current and placed in a moving gas stream, the square of the 
current required to maintain a constant resistance in the wire is a 
linear function of square root of the product of the static pressure 
and the effective cooling velocity of the gas for Mach Nos. less 
than 0.3. 

I^oc(pVo)^ 

If ^ is the angle between the wire axis and the flow direction of 
the cooling gas, then: 

VeS V sin/S 45° ^ /S ^ 135° 

or V sin)3 (l - 1.2(D/L)i cos^yS ) 5° yS < 175° 

.0002 c D/L ^ .001 

where, D/L is the ratio of diameter to length of the wire and V is the 
magnitude of the velocity of the cooling gas. 

Thus if an air stream of known velocity and direction is available, 
a calibration curve may easily be obtained by taking a current reading 
with the hot wire in still air and another with the vdre in the air 
stream. The two points may then be plotted, connected with a straight 
line and used as a calibration curve. Fig. 23. 
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(pVe)i 

Fig. 23. Sample Calibration Curve 
A Wheatstone bridge is normally used to observe the resistance in 
the hot \d.re and maintain it at a constant value. Fig. 24. 




Fig. 24. Anemometer Circuit 

Once a calibration curve is available the hot wire may be used to 
determine velocities of airstreams. It will be noted that for measuring 
velocity with the hot wire, & of 90° should be used. is then taken 
as the velocity because in this orientation the measured velocity is 
very insensitive to directional orientation of the wire. 

If the current is held constant and the hot wire immersed in the 
moving gas stream, the resistance will vary as a function of the product 
of the static pressure and the effective cooling velocity of the gas. 
Fig. 25. 




Fig. 25 . Sample Characteristics Curve 



This fact is useful in determining flow direction. For example, 
at a current of as jQ is varied through 360°, the resistance varies 
as indicated by the heavy lined portion of the above curve while the 
wire is immersed in a constant velocity, constant static pressure stream 
of gas. The minimum value of R would occur when jQ - 90°, 270° and the 
maximum R would occur when /6 - 0°, 1S0°. Another useful observation 
is that the resistance would be the same value for /S > 10° as for 
y6 = 350 °. The galvanometer in the Wheatstone bridge will indicate 
the resistance of the wire. 

It will be noted that the directional sensitivity of the hot wire 
is greatest for near zero for two reasons: 

(1) I dR I is greater for lower 

IdVel 

(2) dVg is maximum near /Q = 0 
d>6 

To measure the flow direction no calibration is required except to 
determine the directional orientation of the hot wire axis with respect 
to some device for measuring angular displacement of the axis such as 
a protractor. 

The hot wire elements are normally mounted on probes across two 
metal posts which are of a material which is readily soldered. Fig. 26. 
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Hi ' 






Hot 

Wire 



Fig. 26. Hot Wire Probe 

This arrangement fixes the wire axis approximately perpendicular 
to the probe axis. A protractor may be fixed to the probe in such a 
manner that, as the probe is rotated about its axis, the angular dis- 
placement from some reference position may be measured. By suitable 
calibration, any given orientation of the \dre axis may be made to 
correspond to a knovm reading on the protractor. 

Using this instrument the flow direction may be determined in the 
plane in which the wire is turned. To accomplish this, the probe is 
inserted into an airstream so that the probe axis is approximately per- 
pendicular to the flow direction and rotated about its own axis. The 
resistance of the hot wire will vary, the maximum value being when the 
wire axis is in alignment with the component of velocity in the plame 
of the investigation. As an example, suppose the maximum resistance 
is observed when the protractor reads approximately 345°. It is 
difficult to tell exactly when the maximum resistance is encountered 
because the galvanometer needle is fluctuating. However, when the pro- 
tractor reads 360°, the galvanometer needle fluctuates about a mark which 
will indicate some slightly smaller resistance. The needle fluctuates 
about this same mark when the protractor reads 331.2°. The protractor 
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reading for a maxlmvun resistance should then be halfway between the 
two readings, or 345 « 6°. 

It may be noted that it is not necessary for the two support posts 
to be the same length, Fig. 26, nor that the wire axis be perpendicular 
to the probe axis. The directional sensitivity suffers when the wire 
axis is not perpendicular to the probe axis, however. 

One possible way of eliminating the effect of perturbation 
velocities would be to use a probe with two hot wires mounted in a cross 
fashion. Fig. 27. The two wires could be placed in different legs 
of the bridge so that effects of perturbations would be cancelling. 

After calibration in a known direction stream it could be used for 
direction determination. This instrument would also have double the 
sensitivity of a single wire. 




Fig. 27. Scheme'- j Using A Two Wire 
Probe For Measuring%low Direction. 
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APPENDIX E 



THE COMPUTER PROGRAM 



APPENDIX E 



THE C0I4PUTER PROGRAM 

A computer program was used to reduce the data and is included 
herein. 

The following symbols were used in the program; 

Inputs 



N 


run nmber 


M 


date 


L 

NN 


month indicator ♦•1 March, -1 April 

■’ V 

position no, where data was taken 


P2 


height of water column in tube to which pressure tap 
at outlet flange is attached - cm. of H 2 O 


PST 


height of water column in tube to which static pressure 
line of flow meter is attached - cm, of H^O 


PATH 


height of water column in tube open to atmosphere 
cm. of H 2 O 


DLP 


difference between total and static pressure in flow 
meter - in. H 2 O 


T2 


total temperature in air at flow meter - °C 


BP 


barometer reading 4, in. Hg 

A 


ALFAl 


one of the angles read at each depth - degrees 


ALFA2 


the other angle read at each depth -.degrees 


TARE 

Outputs 


the negative of flow angle read on the protractor 
at position #2 on the centerline of the inlet duct 


PHI2 


flow coefficient 


PSI2 


pressure coefficient 


Y 


depth 


ANGLE 


flow direction 
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other 



WI 


mass flow rate at beginning of run 


WE 


mass flow rate at end of run 


PSI2I 


pressure coefficient at beginning of nin 


PSI2E 


pressure coefficient at end of run 


PHI2I 


flow coefficient at beginning of run 


PHI2E 


flow coefficient at end of run 


PRI 


pressure ratio at beginning of run 


PRE 


pressure ratio at end of run 


W 


mass flow rate 


PR 


pressure ratio 


ALFA 


intermediate step in averaging process 



Other in 
Subroutine 



P2A 


pressure at discharge flange of pump - PSF 


A2 


2 

discharge area of impeller - ft 


U2 


impeller peripheral speed - ft/sec 


R 


gas constant 


PSTA 


flow meter static pressure - PSF 


TT2 


same temperature as T2 - °R 


BPA 


barometric pressure - PSF 


RHO 


density at flow meter - Ibra/ft^ 


ZETA 


flow rate factor 


T2F 


same temperature etStHZ - °F 


Z 


viscosity of air in flow meter 


VAV 


average velocity in flow meter - ft/sec 


A 


cross-sectional area of discharge pipe - ft^ 


REX 


modified Reynolds number in discharge pipe 


TOEF 


referred flow rate - ft^ 
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For each run, this program requires seven data cards , On the first 
card, the values of N, M, L and NN are entered in four fixed point fields 
of ten. The second and third cards are six floating point fields of 
ten in uhich are entered the values obtained for P2, PST, DLP, T2, PATM, 
and BP» The second card contains those values measured at the beginning 
of the run and the third card, those values measured at the end of the 
run. The next four cards contain eight floating point fields of ten. 

The first of these contains the values of ALFAl(l) through ALFAl(8). 

The second contains the values of ALFA1(9) through ALFAl(15) and TARE. 

The last two cards are the same as cards four and five, except 
substitute ALFA2 for ALFAl. 

The program first enters the subroutine twice to compute flow 
coefficient and pressure coefficient from the data on cards two and 
three. It first changes all pressures to PSFA and temperature to ®R. 

If different manometers and thermometers are used, some of the constants 
may require change. It then does two iterations to find ¥ because of 
the interdependence of W and ZETA. The formula for ZETA in terms of 
REX should not be used for REX greater than 70. Fig. 17. The 
subroutine then computes the pressxire coefficient and the flow 
coefficient. 

The flow coefficients and pressure coefficients at the beginning 
and end of the run are then averaged. The first part of the table is 
then printed. Next the two angles measured for each depth are averaged 
and printed with its corresponding depth. 
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PROGRAM RREROr 

DI MENS I ON Y ( 15) , ALFA! ( 1 5) , ALFA2 ( 15) , ANGLE! 15) 

Y ( 1 ) = 0,. 1 
on 20 J=2, lU 
A = J-1 

20 Y( J) = 0.5' • A 
Y( 15) = 6.9 

25 READ 99, N,M,L,NN 

CALL FLOW (W1,PSI2I,PHI2I,PRI ) 

IF(PRI)7,7,6 

6 CALL FLOW (WE,PSI2.E,PHI2E,PRE) 

W=( WE+WI )/2.0 
PR= ( PRt+PRl ) /2.0 

PS12=(PSI2I+PSI2£)/2.0 , , • 

PHI2=( PHI2I+PHI2E) /2.0 

PRINT no 

PRINT 111, N 
PRINT 1 12 
IF( L)21 ,21 ,22 

21 PRINT 121 ,N,M 
GO TO 23 

22 PRINT 122, N,M 

23 PRINT 123, NN 
PRINT 124, PHI2,PS12 

PRINT 125 . ' 

REAO 300, ALFAl, TARE ■ ' 

READ 300, ALFA2, TARE ' 

DO 30 K=1 , 15 . , • , . 

ALFA= ALFAl (K)+ALFA2(K) + 2.0*TARE 
IF( ALFA-360.0) 1,2,3 

1 IF(ALFA-270.0)4,4,2 

2 ALFA=ALFA + 360.0 
GO TO 4 

3 IF( ALFA-450. 0)5, .5, 4 
5 ALFA = ALFA - 360.0 

4 ANGLE(K)=ALFA/2.0 
30 PRINT 200, Y(K), ANGLE(K) 

GO TO 25 

no FORMAT! IHl/////////) 

111 F0RMAT!38X6H TABLE 13///) 

112 FORMAT! 34X20H PREROTATION TEST ///) 

121 F0RMATM4X8H RUN NO. I3,31X I5,nH APRIL 1964 ///) 

122 FORMAT! 14X8H RUN NO. I3,51X I5,11H MARCH 1964 ///) 

123 F0RMAT!35X13H POSiriON NO. 12///) ■ 

1240F0RMAT! 14X19H FLOW COEFFICIENT = F7.5,. 

13X23H PRESSURE COEFFICIENT = F7.5/// ) , 

125 FORMAT! 32X6H 0£PTH 9X6H ANGLE ///) 

200 FORMAT! 32XF6.2,9XF6.2// ) 

99 F0RMAT!4I10) 

300 F0RMAT!8F10.0) 

7 STOP 
END 

SUBROUTINE FLOW ! W , P S I 2 , PH 1 2, PR ) 

REAO 100 , P2, PST, DLP, T2, PATM, BP 
IF !BP)9,9,B 

8 P2a=! !PATM-P2)/ 34 . 54 U + BP ) *70 . 727 
A2= .1453 
U2= 424.19 
R=53.35 

PSTA=! ! P ATM- PST 1/34.544 +BP)«70.727 
DLPA= DLP *5.2 
TT2 = !T2 + 40. )« 1 .8 +420. 

PPA = BP*70.727 
RHO = PSTA/ ! R* TT2 ) 

PR = P2A/BPA 
ZETA = .7 

VAV = SQRTFI64.343 * ZETA * DLPA/RHO ) 

A = 3. 14159/ 64.0 

W = RHO * A * VAV 
T2F= TT2- 460.0 

Z=.0162 + .000029 *T2F -.00000001 *T2F*T2F 
REX = 2.1053 * W/ Z 
ZETA * .543 + .0093S*REX -. 000 1 65*3 EX*REX • + . 00000 1 *REX* *3 
VAV = S0RTF!64.340*ZETA « DLPA/HHO) 

W= RHO * A * VAV 

WREF = SORTFIR* TT2/32.174) * W /B>A 
PSI2=!P2A-RPA) *32. 1 74 *R»TT2/ ! U2*U2» P2A ) 

PHI2= W«R*TT2/!U2*P2A*A2) 

^ GO TO 11 

9 PR = 0.0 
11 RETURN 

100 FORMAT ! 6F10.0) 

END 
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